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Causality as an emergent macroscopic phenomenon: The Lee- Wick O(N) model 
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In quantum mechanics the deterministic property of classical physics is an emergent phenomenon 
appropriate only on macroscopic scales. Lee and Wick introduced Lorentz invariant quantum the- 
ories where causality is an emergent phenomenon appropriate for macroscopic time scales. In this 
paper we analyze a Lee- Wick version of the O(N) model. We argue that in the large N limit this 
theory has a unitary and Lorentz invariant S matrix and is therefore free of paradoxes in scattering 
experiments. We discuss some of its acausal properties. 



I. INTRODUCTION 



It is interesting to try to understand if one or more 
of the pillars of modern physics may be violated by a 
theory that gives approximately the same experimental 
results as ordinary relativistic quantum field theory for 
experiments that are presently accessible. One such pil- 
lar is causality. Theories that are not causal appear, at 
first glance, to be fraught with paradoxical behavior that 
renders them inconsistent. In the late 1960's Lee and 
Wick [H,[2i proposed an extension of quantum electro- 
dynamics where the Pauli-Villars regulator is treated as 
a finite mass scale. In this theory the Fourier transform of 
the gauge field two-point function has massive "Lee- Wick 
photon" poles with wrong sign residues. It is easy to 
show that this theory is equivalent to a higher derivative 
theory. Naively such theories are unstable and not uni- 
tary. Lee and Wick and Cutkoski, Landshoff, Olive and 
Polkinghorne [3 gave rules (the "LW" and "CLOP pre- 
scriptions") for calculating pcrturbative scattering am- 
plitudes in this higher derivative theory that, for a wide 
class of Feynman diagrams, overcame these difficulties 
yielding Lorentz invariant and unitary scattering ampli- 
tudes. However Lee- Wick electrodynamics is not causal 
for microscopic time scales. Their ideas provide a frame- 
work for studying acausal theories where the acausality is 
only detectable in experiments that can access very high 
energies and/or very short time scales Q. 

In recent papers we extended the work of Lee and Wick 
to non-Abclian gauge theories and argued that they can 
solve the hierarchy puzzle 0, [|| . Even if the ideas of Lee 
and Wick are not relevant for the hierarchy puzzle it is 
worth exploring the physics of acausal theories and exam- 
ining their consistency. Previous work has some limita- 
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tions. While the LW and CLOP prescriptions have been 
shown to give Lorentz invariant scattering amplitudes in 
a large class of Feynman graphs, it is not known whether 
this is true to all orders in perturbation theory. More- 
over, there are serious obstacles to a non-perturbativc 
path integral formulation for the Lee- Wick theory with a 
CLOP prescription ■ Other formulations with prescrip- 
tions different from CLOP's can have a non-perturbativc 
path-integral formulation [g[ , but these have yet to be 
shown to give a Lorentz infant S matrix 0. Perhaps 
there is some subtle obstacle that prevents the construc- 
tion of non-trivial Lee- Wick theories that have a unitary 
and Lorentz invariant 5" matrix to all orders in perturba- 
tion theory. We will argue that this is not the case since 
in large N the Lee- Wick O(N) model provides an exam- 
ple of an acausal theory that has a unitary and Lorentz 
invariant S matrix. 

The leading behavior at large N of the scattering am- 
plitudes in the O(N) scalar model can be calculated fill ]. 
In this paper we consider the Lee- Wick version of this 
theory and argue that at large N the scattering matrix is 
unitary and Lorentz invariant. If the theory has unitary 
time evolution there will be no paradoxes in experiments 
that involve normal scalars in the prepared initial state 
and the observed final state. After all, for any initial 
state there are various possible orthogonal final states 
and the probability for each of them occurring sums to 
unity. This theory has unconventional acausal effects, 
some of which we illustrate with explicit calculations, but 
that does not mean it is inconsistent. 

If a Lee- Wick extension of the standard model is rel- 
evant for the hierarchy puzzle then the acausal effects 
can be studied indirectly in high energy accelerator ex- 
periments through unusual interference effects associated 
with a Lee- Wick resonance [l2[ , like the wrong sign of the 
phase shift of the resonant scattering amplitude. How- 
ever, there is no compelling reason that acausal effects 
should occur at the weak scale; perhaps low energy su- 
persymmctry or a warped extra dimension provide the 
solution to the hierarchy puzzle. There are constraints 
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on the masses of the Lee- Wick resonances from preci- 
sion electroweak physics. These are quite strong because 
integrating out the Lee- Wick resonances gives tree level 
contributions to the oblique parameters S and T [l3j . 

It is possible that acausal effects, of the type we are 
studying in this paper, arise from the extension of the 
standard model to include a quantized theory of grav- 
ity. String theory is an extension of the standard model 
that includes a consistent quantum theory of gravity. At 
the present time it is widely accepted that string the- 
ory is realized in nature. This is because of a lack of 
alternatives and because, even though string theory is 
highly constrained, it has compactifications with enough 
light degrees of freedom to accommodate the known stan- 
dard model physics as well as gravity. However, there is 
no experimental evidence to support the hypothesis that 
string theory is realized in nature. Therefore, even if 
string theory is incompatible with the type of acausality 
we are studying, it seems worth keeping an open mind on 
this issue and contemplating the possibility that acausal 
effects occurring on time scales of order the Planck time 
may arise from the extension of the standard model to 
include quantum gravity p^ |. 

In [l4| it was argued that gravitational radiative cor- 
rections can induce higher derivative terms of the Lee- 
Wick type. In [l5[ a Lee- Wick theory of gravity was 
considered. 

In the auxiliary field formulation of the O(N) model 
at large N the only loop diagrams that enter the calcula- 
tion of scattering amplitudes is the one-loop 1PI auxiliary 
field self-energy. In the Lee- Wick extension of the O(N) 
model this Feynman diagram can be treated using the 
LW and CLOP prescriptions. Hence this theory has a 
unitary and Lorentz invariant S matrix at large N and 
this toy model provides a convenient laboratory to study 
the physics of acausal theories. In this paper we explore 
the Lee- Wick O(N) model. We show by explicit calcu- 
lation that the two particle scattering amplitudes satisfy 
the optical theorem, argue that the S matrix is unitary 
and calculate the acausal behavior that arises in some 
experiments. 

Higher derivative versions of the O(N) model at large 
N have been studied before [lj| and the question of uni- 
tarity is taken up by Liu in Ref. Q. Our investigations 
differ in several respects. The prescription that Liu uses 
differs from CLOP's prescription. In Ref. @ terms of 
order d 6 are added to the Lagrangian, arranged so that 
at tree level there is a complex pair of poles in propaga- 
tors. The imaginary part of these poles is not associated 
with a decay width. In our model only terms of order 
d 4 are added to the Lagrangian so that at tree level the 
two poles in the propagator are at real and positive val- 
ues of p 2 . It is interactions that turn the wrong residue 
pole into a complete "di-pole," and the imaginary part 
is dictated by the physical width. Finally, Liu's proof of 



unitarity is indirect [2p| while we demonstrate unitarity 
by direct calculation. 



II. TIME DEPENDENCE IN A SCATTERING 
EXPERIMENT 

We start our discussion by entertaining the question, 
how does one go about testing causality or looking for 
causality violation in a theory that gives only an S ma- 
trix? One may wonder if this is possible at all since 
formally the S matrix relates states of the infinite past 
to those in the infinite future. Intuitively it is clear that 
this is no impediment: we may prepare two localized 
wave packets in the infinite past to travel toward each 
other and set up detectors to look for the outcome of 
their collision. Moreover, if both the distances traveled 
by the prepared wave-packets to the collision point and 
the distances from this point to the detectors are truly 
macroscopic, then there is an S matrix description of this 
process. Clearly, information on the position and timing 
of first detection of collision products can then give in- 
formation on the causal behavior of the interaction. 

This section formalizes these statements mathemati- 
cally. In theories with normal causal behavior, a reso- 
nant collision that takes place at some space-time point, 
zq, results in the production of outgoing stable parti- 
cles that appear to arise from a second space-time point, 
z' . This second point occurs later in time (t' > to) and 
is typically separated from the collision point by proper 
time of the order of 1/T, the inverse of the width of the 
resonance. The distribution is a decaying exponential. 
This is encapsulated in Eq. (|35|) , in which the separation 
between the two points is w = z' — Zq. 

The situation is quite different for resonant collisions 
through a Lee- Wick resonance. Here the detected parti- 
cles appear to come from z which occurs earlier than zq 
in time (t' < to). The points arc still separated by proper 
time of the order of 1/T, still distributed as an exponen- 
tial that decays away from zq. This is the content of the 
final equation in this chapter, Eq. (|47|) . 

We have been careful to state that the collision "ap- 
pears to" take place at z' . The measurement is made 
a long time after and a long distance from the collision 
region. Within the quantum mechanical S matrix for- 
malism there is no means by which we can investigate, 
nor is there meaning to, the question of precisely where 
or when the collision takes place. This observation is im- 
portant in understanding the interpretation of the results 
in the case where the collision goes through a Lee- Wick 
resonance, where normal causal behavior is violated. 
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Kinematics 



transform 



We prepare from stable particles of mass m an initial 
state consisting of two wave packets traveling towards 
each other from far away. They are initially localized 
about spacetime points y\ and y 2 which we can assume 
are space-like separated, (2/2 — J/i) 2 < and have large 
negative time components (we imagine the interaction 
will take place at around zero time). So we take < 
and \y®\ 3> 1/m where m is the mass of the particles. 
Since they will have to travel a long distance to the in- 
teraction point we also take \yi\ ^ 1/m. 21 

We also want the wave packets to have specific mo- 
menta. That is, their Fourier transforms are localized 
about pi = mvi and P2 = mv 2 . Of course, the momenta 
have to point towards each other so that there is a colli- 
sion. The collision occurs at a point z$ such that 



zo - yi 
n 



v\ and 



zq - V2 

T 2 



I '2 



(1) 



where = yj (zq — yi) 2 is the proper time along the 
world line of the particle from the start point to the in- 
teraction point. 

So we take for the initial state 



IV'in) = / d 4 X! d A x 2 fi(xi - yi)h{x2 - y 2 )4>(xi)4>(x 2 )\0) 

(2) 

with fi(x) concentrated about x = 0, and the Fourier 



fi(h) 



d A xe l 



(3) 



concentrated about ki = pi with pf = m 2 . Here 4>(x) is a 
real scalar field that, when acting on the vacuum, creates 
a stable particle of mass m. 

Similarly, we will set up two detectors for the outgoing 
particles that each record only at a particular point in 
space at a specific time, that is, at spacetime points y[. 
These points can also be taken as space-like separated 
and at late times and large distances, yf 1/m and 
lyll 2> 1/m. And we want to absorb specific momenta, 
p\ = mv[. If the outgoing particles emerge from a point 
z' n then 



and 



!)2 



(4) 



So we take for the final state 



IVw) = / d i x 1 d i x 2 gi{x' 1 -y[)g2{x l 2-y'2)<f>{xi)(j}{x' 2 )\Q) 

(5) 

with giix) concentrated about x = 0, and their Fourier 
transforms concentrated at p[ = mv^. 

Consider now the amplitude for the state \ip- m ) to 
evolve into the state |^ ut), 



(VoutlV'm) = J d 4 xi d 4 x 2 d 4 x\ d i x 2 9*i{x l - y'x)gl{x' 2 - y' 2 )fi{xi - yi)f 2 {x 2 - y 2 ){Q\(j){x' 2 )<p{x l )(j){x 1 )(j){x 2 )\Q). (6) 

Since we have initial points that are space-like separated the order of the fields at x\ and X2 is irrelevant and the 
same goes for the fields at x\ and x' 2 . Also the fields at x\ have later times than those at Xi. So we can replace the 
product of fields above by the time ordered product, which is just the 4-point Green function, 

<0|T[<K4)^K)^i)^2)]|0) = J I] ^^ el(kl ' Xl+k2 ' X2 ~ K< - K - x[ H^) 4 5 {4) (h +k 2 - k[ k' 2 ) 

x TTt2 — \ — -772 — \ — -r (4) (fci,fc 2 ,-fci,-fc 2 ). (7) 
11 k 2 - m 2 + ie k' 2 - m 2 + it K 1 2 ' w 

We have written this in Fourier space in terms of the amputated 4-point function. We will consider cases where 
the amputated 4-point function r^(A;i, k 2 , — k' t , — k' 2 ) is the sum of three terms which depend respectively on the 
Mandelstam variables, s, t and u. At large separation \z$ — Zq| the amplitude (i/'outlV'm) will be negligible except 
when there is a narrow s-channcl resonance in the 4-point function. This is clear if we use customary causal intuition, 
that a narrow resonance can be thought of as a long lived unstable particle produced at zq decaying later at z' , 
but mathematically it is true even when the resonance "decays" at z' before zq- Therefore to examine the leading 
dependence of (VwlV'in) on z ~ z' , for large \zq — z' \, only the term that depends on s = (ki + k 2 ) 2 is important 
and we drop the other pieces. We denote this term by ri^s). 
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Now, multiplying the above by 

1 = / -0^(2nr5(k 1+ k 2 -q) (8) 

we get 

(VwlV'in) = J d A xi d A x 2 d 4 x[ d 4 x' 2 gl(x[ - y[)g* 2 {x' 2 - y' 2 )fi{xi - yi)f 2 (x 2 - y 2 ) 

x e i(k 1 -x 1 +k 2 -x 2 -k' 1 -x' 1 -k' 2 -x' 2 ) TT » * p(4) 

kf — to 2 + ie fe„' 2 — m 2 + ie s 

i 1 1 

= J ^%)<5(g)lfV) (9) 
where we have introduced 

p ^= / n T^w'tfa + fe2 - g ) e <s * , - w /i(fa)^(*a) n fc2 _ ^ 2 + fe (io) 

i i * 

and 

/ n (0( 2 ") 4 ^ fc i + fc 2 - g) e-^ kw *mi)?2(k 2 ) n fc/2 _ v, - te . (id 

i i 1 

The integral in (jlO[) can be broken into two single particle integrals by representing the delta function as an integral, 
Hi) = Jd 4 zJ II fSr j^-^^h{ki)~fM II fc? _ + fe = / ^e^JxWAW (12) 

I 



where 

/lW= /|^/, fe)e »(«,-.»_ F _i 7 _. (13 , 

We now estimate . To this end notice that each com- 
ponent of yi — zq is much larger in magnitude than l/m, 
by assumption. We will compute for z ks zq and come 
back later to see that the integral in the last line of (fl2|) 
has support localized around z — Zq. Rewrite Ii by ex- 
ponentiating the propagator, using the ie of the normal 
particle propagator 

Ii = J dsj ^Mki)e iki < yi - z) e is ^- m2 +^. (14) 



We do first the ki integration by the method of station- 
ary phase. (22| The condition that the phase be stationary 
is 



(ki ■ ( yi -z) + sk 2 ) = 0, (15) 

which implies that 

h = Z -^- (16) 



So we get 



I 
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Now do the s-integration also using stationary phase. The condition is 



which implies that 



and leads to 



d_ 

~d~s 



4s 



s{— m + ie) 



2m 



We now put the pieces together. First, putting the result for F above into F as given in (fT2")) we have 



(18) 



(19) 



(20) 



F ~ / d A ze iz - q 



2 8 7r 3 m 2 (sis 2 ) 3 / 2 



j x ( Z ~VA f 2 ( Z -VA g-imV^-Kil'-im^/Fw) 5 

V 2si / \ 2s 2 ) 



(21) 



r 



Let's investigate this function. Using stationary phase, 
which is justified because we can think of q or m as large, 
or equivalcntly, because we can take H — * 0, we see that 
the stationary phase condition gives z = z* with sat- 
isfying 



q — m 



Hi 



!12 



V( z * - yi ) 2 v(z~^y2y 



= 0. 



(22) 



The second and third terms are just the arguments of the 
functions /j, which, we recall, have localized support. So 
the function _F is non-zero only for 



z* - yi 



= Vl 



and 



Z* - V2 



which implies the support is at 



''2 



2/2 J 



q = mvi + mv 2 = p\ + P2 



(23) 



(24) 



Moreover, we know how to solve for z„ in the region where 
the integral has support: by Eq. ([1]) we see that z* = zq, 
the point where the interaction takes place. We conclude 
then that 



>F(q) 



(25) 



where F is smooth[23j and has support localized at 
q = pi + P2- No te that we are absor bing the p hase 
factor, cxp(-im v /(z - yi) 2 ) exp(-im x /(z - V2) 2 ) = 
exp(— irriTi — irrvr-i), into the definition of F. Similarly, 
we conclude that 

G(q) ~ e- iq <G(q) (26) 
with G smooth and with support localized at q = p[ +p' 2 - 



Using this information on the structure of F and G in 
© we have 



(Vwl^ii 



d 4 q 



-iq<zo-z°)p(q)G(q)TW{q 2 ) , 



(27) 

Eventually we will calculate the w — z' — z dependence 
of (Vwl^in) for large w° in the Lee- Wick 0(N) model. 
Before considering the O(N) model its instructive to con- 
sider some simpler examples. 



B. Standard resonant behavior 

Here we recover the more familiar time dependence as- 
sociated with resonant s-channel exchange in a toy model 
with two real scalar fields <f> and x an d Lagrange density, 



2 M 



(28) 

For simplicity we work at weak coupling g/M <C 1 and 
also take m/M <C 1. The Fourier transform of x's two- 
point function has the form, 



D X (P 2 



ds 



P(s) 



p 2 - s + ie p 2 ~M 2 + iMT 



(29) 

where at order g 2 in perturbation theory the \ width 
is equal to T = g 2 / (32wM). In the narrow resonance 
approximation 



P(s) 



MY 



{s - A/ 2 ) 2 + M 2 r 2 



(30) 
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Since p(s) in Eq. (|30[) is strongly peaked near s = M 2 we 
can extend the s-integration in Eq. ([29]) over the whole 
real s line. Performing this integration using contour 
methods with the value for p(s) in Eq. (|30[) reproduces 
the usual resonance form of the propagator on the far 
right hand side of Eq. ([29]) . 

We calculate the dependence of (?/w |"0in) on w = z' — 
zq for large proper time V w 2 , under the assumption that 



the functions F(q) and G(q) are slowly varying and have 
most of their support around q= p\ +p 2 and q = p[ +p 2l 
respectively. We also assume that the coupling g is small. 
The amputated four point function is, 



(31) 



and so Eq. (|27|) becomes, 



-19 



d A q 

W) 4 



F(q)G(q)e 



1 



q 2 - AP 



iMT 



-g- J ds 
Jo 



d^q 

W) 4 



F(q)G(q)e 



-iqw is(q 2 -M 2 +iMr) 



(32) 

The integration over the components of the momentum q is done using the stationary phase approximation. The 
stationary point is at q = w/ '(2s) and we find that (up to a constant phase), 



(V'outl^ii 



9 



ds l-) F(w/(2s))G(w/(2s))e 



-i{w A / (As)+sM 2 ) e -TMs 



(33) 



Next we perform the s integration using the stationary phase approximation. The stationary point is at s = v 'w 2 / (2M) 
and we arrive at, 



(V'outlV'ii 



g 2 VM 



2(27rvV) 3 / 2 



F(Mw/Vw 2 )G(Mw/Vw 2 )e- 



iM Vw 2 -TVw 2 /2 



(34) 



Since the functions F(q) and G(q) are peaked around q = pi + P2 and q — p[ + p' 2 , respectively, the amplitude is 
appreciable only if p\ +p 2 ~ p[ +p' 2 ~ Mw/Vw 2 . In the center of mass frame, the initial state must be prepared with 
total energy near M and the detectors in the final state are designed to find ordinary particles that are back to back 
with total energy near M. The amplitude is dominated by w° > and w ~ and we can write, in the CM frame, 



(35) 



(VwlV'in) * g(^°) 2( ^ ^3 /2 f(M^/V^)G(Af W /V^)e- a/t "%- r - / 2 



Note that theta function means that the out-going <j> wave-packets appear to emerge from the x decay at a time z' 
that is after the time Zq that the incoming cf> wave packets collide. The factor of exp(— Tw°/2) gives the characteristic 
exponential decay associated with the \ resonance and the factor of (1/w ) 3 ^ 2 arises from the spreading of the x wave 
packet. 



Lee- Wick resonant behavior 



Here we illustrate the acausal behavior of (^out ^in) in 
the simple Lee- Wick toy- model introduced in [jj]. The 
Lagrange density for this theory is, 

C = \d^4> - ^J2(d 2 $) 2 \m 2 ^ 2 - 1# 3 . (36) 

The higher derivative term can be removed by adding a 
field 4> in terms of which the Lagrange density becomes, 



1 



9, 



1 



2 12 

m (p 



f^M 2 4> 2 -^ 3 . (37) 



Next we define <j) = 4>+4> since in terms of <p and <f> the two 
derivative terms are not coupled. The Lagrange density 
now takes the form, 



C = l -d^8^ - \d»m + -M 2 <j> 2 - -m 2 {cj> - 4>f 



1 - r 2 72 1 2/ 



-^9{4>-4>f- (38) 

Provided that M > 2m the mass matrix can be diago- 
nalized by a symplectic transformation 



= cosh 64>' + sinh6»0', 
= sinh6»^' + cosh0(/>' 



(39a) 
(39b) 
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where 

tanh20 = -2m 2 /{M 2 - 2m 2 ). (39c) 
The Lagrange density then takes the form 

c = \d^d^' - l -m' 2 <t>' 2 - \d^d^ + ±Af V 

- ^ 9 (cosh(9-smh6l) 3 (>' -4>'f. (40) 

Defining g' = g(cosh9 — sinh^) 3 and then dropping all the 
primes gives the Lagrange density in a convenient form. 
For simplicity we take m -C M. 

The free field propagator for the normal scalar takes 
the usual form i/(p 2 — m 2 ), however, the free field prop- 
agator for the Lee- Wick field, </>, is —i/{p 2 — M 2 ) which 
differs by an overall minus sign from a conventional scalar 
of mass M. That minus sign means that the propagator 
that one gets from summing Lee- Wick self energy inser- 
tions develops a complex pole at p 2 = M 2 = M 2 + iMT. 
Note that this is has positive imaginary part. Since the 
propagator remains real and regular on a segment of the 
real axis (below the two normal particle cut), the prop- 
agator satisfies Green's reflection principle, (D^(p 2 ))* = 
D^(p 2 *). There is therefore a second pole at p 2 — M* 2 . 
The propagator can be written as the sum of terms with 
poles at M 2 = M 2 +iMT, M* 2 and the two particle cut, 

n 2^ — * — * * f°° j P( s ) 

^ V ] = p 2 - M 2 C + p 2 - M* 2 + n J im , dS p 2 -s + ie- 

(41) 



These poles must not give rise to additional imaginary 
parts in matrix elements since only the normal 4> particle 
is in the spectrum of the theory. 

In the narrow resonance approximation the spectral 
density, p(s), is again given by Eq. ([3H)l . The spectral 
representation for the propagator given in Eq. (|29[) con- 
tains no poles in p 2 but rather has a cut associated with 
the integral over s. However, for a Lee- Wick resonance 
there are poles at p 2 — M 2 and p 2 = M* 2 . In the narrow 
resonance approximation the term in Eq. (|41"j) that has a 
pole at p 2 = M* 2 cancels against the term that contains 
the integral over s (i.e. the cut piece), 

D ^ * p 2 -M 2 -iMT < ™ 

where F ~ g 2 /(32irM). 

We want to calculate the large w° behavior of 
(VWtl^in) that arises from s-channel exchange of the Lee- 
Wick resonance (j> at tree level making the same assump- 
tions that we did in the case where there was s-channcl 
exchange of the ordinary resonance x- I n the case of 
Lee- Wick resonant exchange, 

T^( q 2 ) = -g 2 D^q 2 l (43) 
with given in Eq. (14"2")) . We follow the same steps 
that were used for the \ case. Since the width term in 
the propagator has the opposite sign the phase of the 
exponential proportional to s must be flipped to get con- 
vergence at infinity. Hence, we find that, 



(44) 

The stationary point for the q integration is now at, q = —w/(2s) and we find that (up to a constant phase), 



g 2 r °° /ix2 



(V'outlV'm) - /r) x 2 / ds 



~ (h) H-W(2s))G(-w /(2s))e^ 2 /^ M2 h^ Ms . (45) 
The stationary point for the s integration is in the same place as before, s = \J w 2 / (2M) and we arrive at, 

(VoutlVin) c F(-Mw/V^)G(-Mw/V^y M ^e- r ^/ 2 . (46) 

2(27rv w 

Since in the CM frame the functions F(q) and G(q) are peaked around q° = M and q = the amplitude is dominated 
by w° < and w ~ and we can write 

(Voutl^in) * °(~™ a ) 2{ l ^3/2 H-Mw/Vrf)G(-Mw/Vrf)e- iMw °e Tw °/ 2 . (47) 

The theta function in Eq. (|4T[) means that the out-going <f> wave-packets appear to emerge from the <f> decay at a 
time z' that is before the time zq that the incoming <f> wave packets collide. The factor of cxp(Fu; /2) gives backwards 
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in time exponential decay of the <j> resonance and the factor of (l/|u;°|) 3 / 2 arises from the spreading of the <fi wave 
packet. 



III. REVIEW OF THE O(N) MODEL 

The theory contains N scalar fields 4> a (x), a = 1, . . . N, 
and is invariant under orthogonal transformations of 
them, cj) a (x) — > (j)' a (x) = O ab <j) b . The Lagrange density is 

c = \d,rd»r - - (rrf . (48) 

To get a sensible large N limit of this theory we must 
take Ao ~ 0(1/N). Then the 4>cj> — > 4>(j> scattering ampli- 
tude is 0(1/N), the # -> (jxjxjxj) amplitude is 0(1 /TV 2 ) 
etc.. This is similar to the behavior of QCD in the large 
number of colors limit [l7j . There color singlet mesons 
M with interpolating fields of the form[24[ qq/\/N c have 
MM -> MM scattering amplitudes that are 0(1/N C ) 
and MM — > MMMM scattering amplitudes that are 
C>(l/iV 2 ). However in the case of the O(N) model the 
(jxp — > (jxj) cross section averaged over initial values of the 
O(N) quantum number a and summed over final values is 
(1/ N) while in the large N c limit of QCD MM -> MM 
scattering cross sections are 0(1/N 2 ). 

It is convenient to remove the quartic interaction term 
in Eq. (|48[) by introducing a non dynamical scalar a, and 
make the N dependence explicit by introducing A = iVAo. 
The Lagrangian density takes the form, 

C = \d^ a d^ a - + ^a 2 - \o^^ a . (49) 

One can show that the Lagrange density in Eq. (|49[) is 
equivalent to that in Eq. (|48|) by integrating out the 
scalar a using its equations of motion, a = \(f) a 4> a / (2N) . 
In this formulation the only interaction vertex is between 
a sigma and two </>'s. 

It is straightforward to argue that in this formulation 
of the theory, at large N, the only loop diagrams that 
must be computed is a 1-loop a self energy, So and a 
a tadpole. The physical effects of the tadpole can be 
absorbed into the <j> mass by making the replacement, 
?71q — * to 2 . Using dimensional regularization with MS 

I 



subtraction, 

_ , 2 , N f 1 fm 2 -p 2 x(l~x)-ie\ 

(50) 

where /i is the subtraction point. The sigma propagator 
is, 

DAP 2 ) = ^ + ^MP 2 )^ + ■■■= 77T ,V (n 2 V 

a77 1/\ + U(p) 

(51) 

Writing the scattering matrix as S — 1 + iT, we define 
the scattering amplitude, M, by 

(k[,c; k' 2 , d\T\ki, a; fo, &} 

= (2ir) 4 6 4 (k 1 +k 2 -k' 1 -k , 2 )M(k 1 , a; k 2 , b -> k[,c; k' 2 , d). 

(52) 

Using our expression for the sigma propagator, we find 
that the scattering amplitude is given by 

M(h,a; k 2 ,b -> k[, q k' 2 , d) = ~ ( \ + as(s) + ' ' ') ' 

(53) 

where s, t, u are the usual Mandelstam variables and the 
ellipses denote the two terms similar to the one presented 
that are functions of t and u. Note we have written So = 
TVS to make all the N dependence explicit. 



A. Unitarity of the Two Particle Scattering 
Amplitudes 

With these results in hand, we can explicitly check 
unitarity of two particle scattering in the O(N) model to 
leading order in 1/iV and to all orders in A. Unitarity 
of the S matrix, i.e., S^S = 1, is equivalent to i(T^ — 
T) = T^T . Taking the two particle matrix clement of 
this equation gives 



i (M(k[,c; k' 2 ,d— > ki,a; k 2 ,b)* - M(h,a; k 2l b-^ k[,c; k 2 ,d)) = y^M(fci,q; k 2 ,b -> ip)M* (k' v c; k' 2 ,d^ ip). (54) 

■4' 

To simplify this equation, note firstly that at leading order in N , we may restrict the summation above to two particle 
states. Next, we use the fact that the theory is invariant under the combined time reversal times parity discrete 
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symmetry, so that M.(a — ► b)* = A4(b — > a). Thus, the requirement of unitarity becomes 

d 3 qi _ d 3 q 2 
< bef I (I-k^IE. (1 

ej 



2ImM(h,a;k 2 ,b^k[,c;k' 2 ,d) = Y,<<ef J j££^ 1 {2 tf2E 2 {21 



x M{k ll a\k 2l b -> q±, e; q 2 , f)M*(k{, c; k' 2 , d — > qi,e;q 2 ,f). (55) 

where the identical particles factor t e / is equal to 1/2 if e = / and is unity otherwise. It is now trivial to check 
unitarity. First, notice that we can express the one- loop correction £(s) as 



nv 2 ) = 



i 



32tt 2 



/ dx log |m 2 — p 2 x(l — x) 
Jo 



/ Am 2 
inJl- -^-d{p 2 -Am 2 ) 



Therefore, since both t and u are negative, the imaginary part on the left hand side of Eq. 

A 2 



is given by 



2ImA^(fcia; k 2 b — > k[c; k' 2 d) 



1 — ^rSabScd 



16ttN |l + AS(s)r 



(56) 



(57) 



On the other hand, the sum over e, / on the right hand side of Eq. (|55p is enhanced by one power of TV when the 
scattering is in the s channel. Thus, to leading order, the sum is given by 



dhl 



d\ 2 



(2tt) 4 (5 4 ((7i + q 2 - pi - p 2 )M(ki,a; k 2 ,b -> q lt e; q 2 , f)M*(k[,c; k' 2 , d -> ?i, e; q 2 , /) 



d 3 fc 2 



2 Qb cd ./ (2tt) 3 2£;i (2^) 3 2£ 2 



(2 7 r) 4 S 4 (k 1 + k 2 -p) 



A 



iVl + AE(s) 



A 2 



16nN 



SabS, 



cd 



Am 2 



1 



l + AS(s) 



(58) 



where s = p 2 . Consequently, we see that the S'-matrix of the theory is unitary to leading order in TV on the two 
particle subspace of the Hubert space. Notice that this argument was sensitive only to the imaginary part of £(p 2 ). 
We will see that in the Lee- Wick case, the real part of the one-loop correction is changed but the imaginary part (for 
two particle final states) is the same. Since the only nontrivial imaginary part is associated with the a propagator it 
should be clear that unitarity also holds for the higher parti cle parts of the Hilbert space . 



Time dependence of Two Particle Scattering 
Processes 



In preparation for our discussion of acausal processes in 
the Lee- Wick O(N) model, we review some aspects of the 
time dependence of the two particle scattering amplitude 
in the normal O(N) model. To simplify the discussion 
we will assume A C 1 and work to one-loop order in 
perturbation theory. This approximation allows us to do 
explicit computations while retaining the salient features 
of the causal structure of the time dependent amplitude. 
Using the results of section !!! A[ the transition amplitude 
(VwlV'in) is given by Eq. (g7]): 



where w = z — zq. The four point function (q 2 ) 
can be deduced from Eq. (|53p by expanding in the small 
parameter A. We ignore the tree-level amplitude as 
it leads to trivial time dependence of the amplitude, 
exp(— i(p\ + p 2 ) ■ w) times a function localized about 
w = 0. The one-loop four point function describing 
(a, a) — > (&, b) scattering is given by 



rfV) 



-iX 2 



nl , , ( m 2 — q 2 x{\ — x) — it 
ax log 



d 4 q 



i ^F(q)G(q)T^(q 2 ) 



(59) Thus, the transition amplitude is 



(60) 



(V'outlV'in) = 



-lX 2 

32tt 2 N 



dx 



d 4 q 
(2nY 



F(q)G(q)e^ 



{-tw° dq°) (' 



lq w log(m 2 - q 2 x(l -x)- ie), (61) 



10 



where we have introduced a derivative with respect to q°. Integrating by parts, this derivative acts on F^-q), G(q) and 
the logarithm. Since the functions F (q) and G(q) are slowly varying, we will only keep the term where the derivative 
acts on the logarithm. Therefore the amplitude can be written as 

(Vwl^in) * ~ 2 X l f dx I 7 ^F(q)G(q)e- i " w - 2 2 f (62) 

32tt z Nw v Jq J (27r) 4 q z — m z (x) + te 



where m(x) = m/ \/x(\ — x). Introducing an integration over a variable s to write the propagator as a phase gives, 

(^tl^) * I dx I ds J ^(^(^--e^-^^V. (63) 

As before, we use the stationary phase approximation to evaluate the various integrations. The stationary point 
for the integrations over the components of q is located at q = w /(2s) and performing these integrations gives (up to 
a constant phase) 

\ 2 i pi poo / -i \ 3 

(VwhM^ 16n 2 N {2n) 2 I dx I ds {y s ) e- i ^ 2 /^ +m2 ^ s - i ^F(w/(2s))G(w/(2s)). (64) 



Next the s integration is performed using stationary phase. The stationary point is at s = v « ,2 /(2m(.r)) and we find 
that 

(VoutlV-in) * 3 * 2N (^f) 5/ jyx(m(x)) 3 / 2 F(m(x)w/V^)G(m(x)w/V^)e- im ^^. (65) 

Finally, we have to do the x integral. We use the method of stationary phase once again. The stationary point is at 
x = 1/2 and the transition amplitude is 

(VwlVm) = -g=- e - 2imV ^F (2mw/Vw*)G(2mw/Vrf). (66) 

Recall the functions F(q) and G(q) have support at q = pi + p 2 and q = p[ + p' 2 , respectively. The amplitude is 
appreciable only p\ +P2 ~ p'i +p' 2 ~ 2mw/y/ w 2 . Since the energy is positive, and choosing the CM frame, the above 
can be rewritten as 

It is worth comparing this expression with the transition amplitude in the case where the scattering is mediated by 
a resonance, Eq. (|35[) . In that case there an exponential decay due to the width of the resonance, which is absent in 
Eq. (|6"T|) because the mediators are stable. In the tree-level case the scattering is mediated by one particle; its wave 
packet spreads out like (w )'. In the loop case the presence of two wave packets leads to a powerdaw fall off of the 
amplitude as (w ) 3 . The ^-function in Eq. (|6T[) indicates that the decay particles appear at times after the collision. 



IV. THE LEE- WICK O(N) MODEL 

Let us now move on to study the Lee-Wick O(N) 
model. We begin by discussing the Lagrangian of the 
model before moving on to examine the loop structure. 
Once the loop structure is understood at leading order 
in the 1/7V expansion, we will compute the two particle 
scattering to leading order in 1/N and to all order in A. 
We will use these results to demonstrate unitarity of the 



theory. Finally, we will explicitly compute the time de- 
pendence of one-loop scattering processes, demonstrating 
aspects of the acausality of the model. 
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A. The Lagrangian 

The theory is the usual O(N) model, augmented with 
a higher derivative term. The Lagrangian is given by 

(68) 

where, as in the normal O(N) model, the fields cf> a are 
scalar fields in the fundamental representation of the 
group O(N). We can remove the higher derivative term 
from the Lagrangian by introducing N scalar fields (j) a . 
Then an equivalent Lagrange density is 



the Lagrangian may be written as 

- \d^ a d^ a + \m 2 ff - \a (f - r) (f - 4> a ) ■ 

(72) 

In this formulation of the theory, it is straightforward to 
establish a power-counting argument which shows that at 
leading order in 1/N, the only relevant graphs are the a 
tadpole and the 1PI a self-energy Y,o(q 2 ). We can absorb 
the effects of the a tadpole by replacing TOq by m 2 . In 
the following, we will treat to as a small parameter. It 
remains to compute the self-energy, So( ( ? 2 ) of a to leading 
order in 1/N. 



+ -M 2 d> a d> a - — ' 

2 Y Y 8N 



(69) 



We may diagonalize the derivative terms by defining <f> a = 
cj) a + (p a and performing an integration by parts. The 
Lagrangian becomes 



2 ^ 

\m 2 



8N 



2" MS 

2 



(70) 



This Lagrangian has a simple interpretation. There are 
N normal scalar fields 4> a and N Lee- Wick scalar fields 
(j) a ; these fields have quartic interactions. Note that the 
Lee- Wick scalars (f> a can decay to three 4> a quanta pro- 
vided that (p a is heavy enough. We will assume this 
decay channel is open so that the width T of the Lee- 
Wick scalars is non-zero. There is mass mixing between 
the normal and Lee- Wick scalars which can be removed 
by a symplectic transformation on the fields or treated 
as a perturbation when the ordinary scalars are very 
light compared with the Lee- Wick scalars. Neglecting 
the mass of the ordinary scalars, for large N and small 
coupling A, 



X 2 M 
2 l0 n 3 N' 



(71) 



To obtain a unitary S matrix it is crucial that the <f> 
propagator have poles at complex energy, as in (|41|) . For 
this reason, even though the width is order 1/N it will 
be important to retain it in intermediate steps in our 
calculations. 

As was the case in the normal O(N) model, it is conve- 
nient to introduce a non-dynamical scalar field a so that 



B. Loops 

Before we embark on the computation of the self- 
energy, let us pause for a moment to consider the prop- 
erties of the Lee- Wick resonances. We are familiar with 
the properties of ordinary resonances in quantum field 
theory. One familiar fact is the importance of resuming 
the width of a resonance in order to avoid the appearance 
of spurious poles in Feynman graphs. This resummation 
changes the analytic structure of the theory in a man- 
ner consistent with the non-perturbative information of 
the Lehmann representation. Similarly, resummation of 
the width of Lee- Wick resonances changes the analytic 
structure of the theory in crucial ways. 

At tree level, the Lee- Wick propagator is 



b(j?) 



M 2 ' 



(73) 



At loop level, the particle develops a width. Just as in 
the example we discussed earlier, in Sec lII CI the loop- 
corrected Lee- Wick propagator is given by 



D(p 2 ) 



P 



M* 



ds- 



P(s) 



s + le 
(74) 

where M 2 = M 2 + iMT. We shall use this form of the 
propagator to compute the a self-energy even though the 
corrections due to the width are formally subdominant in 
the 1/N expansion. The subdominant corrections modify 
the analytic structure and it is this modification that 
allows the theory to be unitary. 

The poles present in the Lee- Wick propagator are in 
unusual locations in the complex p 2 plane, so we must 
take care to define the contour of integration in Feynman 
graphs appropriately. We must understand how to define 
expressions such as 



I = 



(2tt) 4 (p + q) 2 - M 2 p 2 



Ml 



(75) 
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where Mi and Mi may be complex masses, either in the 
upper or lower half plane of the Feynman integration. 
Let us consider the p° integral. The integrand has four 
poles. Two of the poles are located at 



V 



Ml 



(76) 



The location of the other two poles depends on the value 
of the external four momentum q. For time-like q we can 
go to a frame where q = and these two poles are located 
at, 



p° = -q° ± 



P 2 



Ml 



The contour Lee and Wick suggested is such that, once 
the Green's functions are computed by Fourier transform 
from momentum space to space-time, there is no expo- 
nential growth in time, and can be described as follows. 
Consider the position of the poles as a function of the 
coupling A present in the theory. At A = the widths 
vanish, so M\ and M2 are real masses. Then the con- 
tour is defined to be the usual Feynman contour. As A 
increases away from zero the Lee- Wick particles become 
unstable, the poles on the real line become complex pairs 
of poles that move away from the real axis. The Lee- Wick 
prescription is to deform the contour, as A increases from 
zero, so that the complex poles do not cross the contour; a 
pole which was initially below the contour remains below 
the contour, for example. If the external momentum is in 



Im p" 



Rep" 



FIG. 1: Contour given by the Lee- Wick prescription for in- 
tegration in the complex p° plane. The crosses denote the 
poles at p° = ± \/p 2 + Mi and at p° 
circles those at p° = 



+ M* 2 and the 
The heavy line 



denotes the cuts on the real axis starting at ±3m. The con- 
tour of integration is deformed as the interactions are turned 
on and the LW poles move into the complex plane so that the 
complex poles do not cross the contour. 



\q°\ < I Mi +M2I. As the momentum is varied (for fixed, 
non-zero A) poles may cross a contour. However, the in- 
tegral can still be defined by deforming the contour so as 
to avoid the pole. This leads to a well-defined contour 
unless poles pinch it. The pinching occurs when a pole 
in Eq. (|76|) coincides with one in Eq. (|77|) , and signals 
the presence of a singularity, usually a branch cut, in the 
integral / (as a function of q). An additional prescription 
is required to define the integral in this case. 

We need a prescription only when the new singularity 
occurs for real valued energy q°. This may occur if one 
propagator carries mass M 2 while the other one has mass 
A/ n * 2 — that is, we could have M 2 = Af n * 2 in Eq. ([Toll . 



^ ' Then it is easy to see that when q° satisfies the equation 



(q°f = 2{f + ReMi 2 ) + 2\f + M 2 



(78) 



two of the poles in Eq. (|76[) and Eq. ([77]) overlap and 
the contour is pinched. The CLOP prescription is as fol- 
lows: Define the Feynman integral by taking the masses 
M 2 and M± 2 to be unrelated complex mass parameters 
so that the poles do not overlap. At the end of the cal- 
culation impose the condition that M| is the complex 
conjugate of M 2 . With this prescription, the self-energy 
So is unambiguously defined and Lorentz invariant, and 
may be computed using standard methods. In particular, 
the contour we have chosen allows us to Wick-rotate the 
integral in Eq. (|75[) and in all Feynman integrals we will 
encounter. In the following, we will compute the integrals 
in dimensional regularization and discard the divergent 
pieces that are proportional to l/(d — 4) in addition to 
the finite pieces involving the logarithm of An and Eulcrs 
constant. In the Lee- Wick O(N) model they cancel, since 
the theory is finite by naive power counting in the higher 
derivative formulation. 123 



C. Computation of the Self-Energy 

We now embark on the explicit computation of the self- 
energy. We define £ = N'Sq as before. There are various 
graphs contributing. The graph involving only the nor- 
mal particles reproduces the self-energy in the normal 
O(N) model, Eq. (f56|) . To simplify the notation, we de- 
fine a function 



F(M 2 ,M 2 ,q 2 



d d p 



(2ir) d (p+q) 2 -M'f 



p 2 - M 2 ' 



(79) 



the unphysical region this prescription is unambiguous. 
For example, for the integral in (|75p . one can start with 



Discarding the divergent piece that is proportional to 
l/(d— 4) and the finite constant pieces involving the log- 
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arith.ni of 4tt and Euler's constant we find 

F{MlMW) = -^ 
x f dx log (xM 2 + (1 - x)M$ - x(l ~ x)q 2 ) . (80) 

It is easy to verify that the only candidate singularity of 
F(M 2 , M|, q 2 ), as a function of the complex variable q 2 , 
is a branch cut with branch point at q 2 = (Mi + M2) 2 . 
The tr self-energy can then be expressed in terms of sums 
of these functions evaluated at various arguments. Thus, 
the contribution, Si of the normal particles to the total 
self-energy, £, is 

^ 1 (q 2 ) = F{m 2 ,m 2 ,q 2 ). (81) 



Since one of our main goals is to understand unitarity of 
the theory, we will focus on understanding any possible 
imaginary parts of £. 



Next we consider graphs with both propagators being 
of the Lee- Wick field, (|?4")) . It is convenient to consider 
the contributions of the pole parts of graphs involving 
the Lcc-Wick particles separately from contributions in- 
volving the spectral density p(s). Firstly, consider the 
terms involving only the Lee- Wick poles. Following the 
CLOP prescription we use different complex masses for 
the two propagators, with M 2 = M 2 + iS and M 2 = AI 2 . 
The loop integral is 



S 2 (g 2 



d d P 



1 



M 2 - id 



- Ml 



id 



_(p + q) 2 -M 2 (p + q) 2 - 
F(M 2 + id, M 2 , q 2 ) + F(M* 2 - id, M* 2 , q 2 ) + F(M 2 + id, M* 2 , q 2 ) + F(M* 2 



Ml 



id, Ml, q 2 ) 



(82) 



It is easy to see that £2 is continuous across the real line. The CLOP prescription has effectively moved the two 
branch points that would have occurred at q 2 = (M c + M*) 2 away from the real axis, by an amount of order 6, to 
y M 2 + iS + M* and y 1 M* 2 — iS + M c . The remaining two terms appearing in the self-energy have complex branch 
points even for (5 = 0. Hence the discontinuity across the real line vanishes, and this persists in the limit that 6 
goes to zero. An explicit example may clarify this. The expression F(M% + id, M* 2 ) + F{M* 2 — id, M 2 ) contains the 
dangerous terms appearing in the Feynman integral in which poles on opposite sides of the contour may pinch when 
5 = (and q 2 is real.) But this expression is explicitly real on the real axis and analytic in a band of width ~ <5 
containing the whole real axis: 

F{M 2 C + id, M* 2 , q 2 ) + F(M* 2 - id, M 2 C , q 2 ) 

dx [log(a;(M c 2 + id) + (1 - x)M* 2 - x{\ - x)q 2 ) + \og(x(M* 2 + id) + (1 - x)M 2 c - x{\ - x)q 2 )] 



32tt 2 



cfalog |z(M 2 + id) + (1 - x)Mf - x(l - x)q 2 



(83) 



Since for real valued q 2 the imaginary part of this vanishes (there is no need to define this as a discontinuity) 
independent of <5 ^ 0. In the limit as (5 — > the imaginary part remains zero. In the remainder of this section, we will 
omit the parameter <5 to simplify the equations. 

In the next section we will use the self-energy to verify the unitarity of the S matrix in this theory. It will be 
useful to write the result for the self-energy concisely While the width is of order l/N, its presence is crucial in 
demonstrating that £2 is real. But once we established that the CLOP defined £2 is real we can neglect the width 
and give a simple expression for the self-energy: 



1 f 1 

E 2(<7 2 ) =-T7— 2 \ dx\og\M 2 -x{\-x)q 2 \ 2 . 



(84) 



Next, we compute the Feynman integrals involving the Lee- Wick pole and the normal particle. We find that the 
self-energy in this case is given by 



£ 3 (s 2 ) = \ J 



d d p 
(2n) d 



1 



1 



1 



m 2 \{p + q) 2 - Ml (p + qf 



Ml 



1 



(p + q? 



1 



1 



1 jy^ 1 



Ml p 2 



Ml 



-2F(m 2 , M 2 c ,q 2 ) - 2F(m 2 , M* 2 ,q 2 



(85) 
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The branch points are both off the reai axis, at q 2 = (M c + m) 2 and (M* + to) 2 . As above, since S3 is reai we can 
neglect the width and write, concisely: 

1 r 1 

£ 3 (<Z 2 ) = T7TT / dxlog\xM 2 + (l-x)m 2 ~x(l~x)q 2 \ 2 . (86) 

Now we turn to terms involving the spectral density p(s). Terms in the self-energy involving products of the 
Lee- Wick poles and the spectral density give 



£4 = - / dsp(s) I 

K J9m2 J 



d d P 1 / 1 1 



(27r) d p 2 -s + ie \(p + q) 2 - M* 2 (p + q) 2 - M 2 , 
dsp(s)[F( S ,Mlq 2 )+F(s,Mf,q 2 )] . (87) 



The sum of F functions is similar to that appearing in Eq. (|S5|) : evidently this is also real and so the integral against 
p is real. In the narrow resonance approximation, we find that self-energy due to these terms is 



1 r 1 

S4= 16^2 7 dxl °s\M 2 -x(l-x)q 2 \ 2 , 



(88) 



Finally, there are terms involving the spectral density p and the normal pole, and involving a double integral over 
two powers of p. These terms do lead to an imaginary part, describing real scattering from two particle states into four 
or six particle states. It is important to understand that the imaginary parts arising from these expressions involve 
final states containing only normal particles. In the narrow resonance approximation, we find that these terms lead 
to a contribution to the self-energy given by 



1 C i-K / AM 2 

— 2 I |M= - x(l - + —^1 - - 4M*) 

iclogliJU 3 + (1 - x)m 2 - x(l - x)q*\ + ^(1- {m + U)2 ) 8(,» - (771 + M) 2 ). (89) 

In total, we find an explicit expression for the self-energy which is simple when we treat the width T to be negligible 
compared to the mass M: 



nf) = - 



32tt 2 



1 \x(l-x)q 2 \\M 2 -x(l-x)q 2 \ 2 
, dxl ° g \xM 2 - x(l - x)q 2 \ 2 l7T9{q ) 

I AM 2 . . / M 2 \ 

, (90) 



_ i7r J l _ ^0(q 2 - AM 2 ) - 2it: (l - 6(q 2 - M 2 ) 

where we have neglected the normal mass m. Notice that the width T of the Lee Wick resonances does not appear in 
this result. It was important for defining the contour for the loop integration but not in the final form of the answer. 
It will also play a role in our understanding of the unitarity of the S matrix. 



D. Unitarity 



Unitarity of the S matrix is equivalent to requiring i(T* — T) = T^T . We consider two particle matrix elements of 
the right and left hand side of this equation and verify their equality to leading order in l/N. For convenience, we 
restate the requirement of unitarity for the amplitude describing scattering of a two particle state into a two particle 
state: 

i (M(k[,c; k' 2 ,d^ fci,a; k 2 ,b)* - M(h,a;k 2 ,b -> k[,c; k' 2 ,d)) = y^M(fc 1 ,a;fc 2 ,b-> ip)M*(k[,c] k? 2 ,d-np), (91) 

■4, 



15 



where M(ki, a; k 2 , b — ► k[,c]k' 2 ,d) is the amplitude for the two particle scattering. Previously, in our discussion of 
unitarity of scattering in the normal 0{N) model, we argued that the only allowed final state at leading order is 
the two particle final state. However, the situation is different in the Lee- Wick 0(N) model. At leading order, two, 
four and six particle final states are accessible. Intuitively, this is because we can create Lee- Wick resonances which 
subsequently decay into three normal particles with unit probability. This is the reason it was necessary to retain 
the width of the Lee- Wick resonances, even if it is subleading in 1/N; after all, for a non-zero width, however small, 
given enough time the unstable "particle" will decay. Therefore we will have to include these additional final states 
in the sum of the right hand side of Eq. (j9"Tj) . In Eq. (pTTj) the initial and final states and the intermediate states ip 
only involve the stable ordinary particles. Even though the propagator for the Lee- Wick resonances contains poles at 
p 2 = M 2 and p 2 = M* 2 these particles are not considered to be in the spectrum of the theory. 

The two particle scattering amplitude in the Lee- Wick theory is given in terms of the self-energy by the same 
expression as in the normal O(N) model: 
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(92) 



where s = (ki + k^) 2 as usual, and the dots indicate the t and u channel terms in addition to higher order terms in 
1/N. Since t and u are negative quantities for physical scattering, we find that the left hand side of the unitarity 
relation Eq. ([91]) is 
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(93) 



For simplicity, and without loss of generality, here and below we neglect the normal mass m. Now we must compute 
the right hand side of the unitarity relation. This is straightforward when the state \4>) in Eq. (|9"Tj) is a two particle 
state; in that case the sum becomes an integral over the two body phase space of the normal particles and the 
amplitude is simply the two- two scattering amplitude of Eq. (|92p . Since we are neglecting the mass of 4> a the sum 
becomes 
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We see that this part of the complete sum over accessible final states \tp) reproduces the first term in parentheses 
on the right hand side of Eq. . The other two terms arise from four and six particle final states. Consider first 
the process with a four particle final state. To leading order in N, the two-four particle scattering must contain an 
intermediate <f> which then decays; the Feynman graph is shown in Fig. O 





FIG. 2: Feynman diagram for the two to four scattering amplitude proceeding through the decay of a Lee- Wick particle. The 
solid line denotes a "normal" particle, the zig-zag line a "Lee- Wick" particle, and the dashed line with the shaded blob denotes 
the dressed a auxiliary field propagator. 



Since we are summing over all final states in the decay of the Lee- Wick resonance, the computation reduces to 
computing the amplitude to make the intermediate state containing a Lee- Wick resonance and one normal particle, 
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and integrating over two particle phase space. This is familiar from the theory of ordinary resonances. In more 
detail, the four particle phase space integral can be organized into an integral over the three body phase space for the 
Lee- Wick decay and an integral over the two particle phase space of the Lee- Wick and fourth normal particle times 
an energy-momentum conserving delta function. The result of the three particle phase space integral is the width of 
the Lee- Wick particle; however, in the region of the two body phase space where the Lee- Wick is nearly on-shcll a 
factor of 1/r appears from the Lee- Wick propagator. Thus, the whole computation reduces to a simple integral over 
the two body phase space. The remaining parts of the integrand can then be interpreted as the amplitude to create 
the intermediate Lee- Wick particle plus one of the final state normal particles. 

The amplitude to create the intermediate state is given up to sign by Eq. (|9"!2)) . Thus, we can explicitly perform the 
sum over four body phase space to find 
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Similarly, the six particle phase space integral becomes an integral over the two particle phase space of two intermediate 
Lee- Wicks. Near the region where the Lee- Wicks are on-shell, there is an enhancement by 1/L 2 . The result is 
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Since the sum of Eqs. (|M|) . (|93|) and (|56"|) equals Eq. (|9"5|) we have verified, for two particle matrix elements, that 
i(Tt —T) = T^T. 

It is easy to extend this argument to show unitarity for any matrix element. To leading order in 1 /N any amplitude is 
given by a sum of skeleton diagrams with the propagators, including the a propagator, replaced by the full propagators. 
In the absence of a Kallen-Lehman decomposition, we cannot proceed with the usual cutting rules to show unitarity. 
For example, it is not obvious, if at all, how to set up a "largest time equation [l8j]. However, one can still analyze 
individual graphs by cutting the diagrams. A cut through a tr propagator is handled using the results for the 2 — > 2 
amplitude demonstrated above. Cuts through normal particle propagators never produce an imaginary part: they 
are just as in the standard analysis and, since we only have skeleton graphs, these propagators are never on-shcll. 
Finally there are "cuts" through the Lee- Wick propagators. These just correspond to taking the imaginary part of 
D(p 2 ) in (|74[). The imaginary part of the sum of complex poles vanishes. We are left with the imaginary part of 
the integral over the spectral function p(s). This has precisely the structure that a normal resonance in the standard 
unitarity analysis has, so it leads to the correct unitarity relation. In particular, it corresponds to a sum, in T^T, over 
intermediate three normal particle states. 

Since the S matrix provides a one-to-one map from the past to the future in scattering experiments, the existence of 
a well-defined S matrix is enough to show that there can be no paradoxes in these scattering processes. Nevertheless, 
the theory is acausal as we shall now explore. 



E. Time Dependence: Acausality 

To study the time dependence of scattering in the Lee- 
Wick theory, we will work to one-loop order in pertur- 
bation theory. The graph containing normal particles 
reproduces the transition amplitude of the normal O(N) 
model, shown in Eq. (|67[) . Our main focus is on the 
acausal behavior associated with poles in the upper half 
plane. All the acausality decays exponentially with time 
except for the case where in the loop one of the poles 
is at M 2 and the other at M* 2 . Then for real incoming 
momentum one can create an on shell configuration with 
two Lcc-Wick resonances and this leads to acausal be- 
havior that falls off with a power of time. In this section 
we calculate this power law acausal behavior. The part 



of the four point function with Lee- Wick poles at M 2 and 
M* 2 is 

r W) = TTrft? f dx\ g{M 2 -i{l-2x)MV-x{l-x)q 2 ). 
1Q-k 2 N J 

(97) 

Since the sign of the Mr term in the logarithm changes 
sign over the region of integration, it is convenient to 
break the integral into two terms as 

2 - 

If V) = Y^JJ dx [\og{M 2 - i(l - 2x)MT 

- x(l - x)q 2 ) + log(M 2 + i(l - 2x)MT -x(l - x)q 2 )] . 

(98) 
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We make the same assumptions as before. In particular as 
the functions F(q) and G(q) are taken to be slowly vary- 
ing and to have support around q° = 2AI and q = 0. We (VwlV'in) = (V>out|V'in)+ + (VWi#in)- (99) 
find it convenient to decompose the transition amplitude 



where 

\ 2 rh r . . O o 

(100) 
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Here, M(x) = M/ ' x{\ — x). We now put the denominator of the propagator into an exponential by introducing an 
integration over a proper time variable s, 

±i\ 2 f°° f d 4 a 

(V>out|V>in>± - 8n 2 Nw o J o dx J o ds J J^4 F (l)G(q)q° e^p (-iq ■ w ± is{q 2 - M{xf) - s(l - 2x)MT/(x(l - a;))) . 

(101) 

It is now straightforward to successively do the q, s and x integrations using the stationary phase approximation. The 
stationary points are at q = ±w/(2s), s — v 'w 2 / (2M(x)) and x = 1/2 and we find that (up to an overall constant 
phase) 

(Voutl^i ^zg e^ 2M ^F(±2Mw/V^)G(±2Mw/V^). (102) 

32VW 2 tt 3 N 

Given where the functions F and G have support we can rewrite this as, 

(VwlVWi - d -^^^^ aMw °H^Mwl^)G{±2Mw/Vw^). (103) 

We have checked by explicit calculation that the other one loop contributions arc exponentially suppressed in |u>°| and 
so for very large the power law term that falls off as l/|iu°| 3 , displayed above, dominates the acausality in the 
one-loop contribution to (Vwl^in)- Note that Eq. (|103[) has a very different behavior than one would expect based 
on the example of single Lee- Wick resonant exchange that we discussed earlier. It is not exponentially suppressed for 
large times and contains both acausal and causal pieces. 



V. CONCLUDING REMARKS 

We have studied the Lee- Wick 0(N) model and argued 
that the prescription of Lee and Wick and Cutkowski et. 
al. yields an S-matrix for this theory that is unitary 
and Lorentz invariant in large N. This suggests that, 
even though the theory is not causal, there will not be 
paradoxical behavior in scattering experiments. 

In this model we demonstrated, by explicit calculation, 
some of the acausal behavior in two-two scattering of the 
ordinary scalars that arises from virtual "Lee- Wick parti- 
cles." The Lee- Wick 0(N) model presents a playground 
to examine the consistency of theories where causality 
emerges only for long enough times and low enough en- 



ergies. There are other theories that are worth exploring 
for this purpose. For example, there are two dimensional 
models that can be solved exactly and it would be inter- 
esting to see if Lee- Wick versions of some of these theories 
are also soluble and, if so, explore their properties. 

Acknowledgments 

We thank David Politzer for important contributions 
to this work. The work of BG, DO and MBW was sup- 
ported in part by the US Department of Energy under 
contracts DE-FG03-97ER40546, DE-FG02-90ER40542 
and DE-FG03-92ER40701, respectively. 



[1] T. D. Lee and G. C. Wick, Nucl. Phys. B 9, (1969) 209. [2] T. D. Lee and G. C. Wick, Phys. Rev. D2, (1970) 1033. 



18 



[3] R. E. Cutkowsky, P. V. Landshoff, D. I. Olive and 
J. C. Polkinghorne, Nucl. Phys. B 12 (1969) 281. 

[4] S. Coleman, "Acausality," In Erice 1969, Ettore Majo- 
rana School On Subnuclear Phenomena, A. Zichichi, Ed- 
itor, Academic Press, New York, 1970, p. 282. 

[5] B. Grinstein, D. O' Connell and M. B . Wise, Phys. Rev. 
D 77 (2008) 025012 |arXiv:0704. 18451 [hep-ph]]. 

[6] B. Grinstein, D. Q' Connell and M. B . Wise, Phys. Rev. 
D 77 (2008) 065010 [arXiv: 0710. 55281 [hep-ph]]; J. R. Es- 
pinosa, B. Grinstein, D. O'Connell and M. B. Wise, 
Phys. Rev. D 77 (2008) 085002 |arXiv:0705. 11881 [hep- 
ph]]; B. Grinstein and D. O'Connell, [arXiv:08 01~I034l 
[hep-ph]]; T. R. Dulaney and M. B. Wise, Phys. Lett. B 
658 (2008) 230 |arXiv:0708.0"567l [hep-ph]]. 

[7] D. G. Boulware and D. J. Gross, Nucl. Phys. B 233 
(1984) 1. 

[8] A. van Tonder, Int. J Mod. Phys. A 22, 2563 (2007) 

[arXiv:hep-th/0610185 . 
[9] C. Liu, "Strongly Interacting Higgs Sector Without 

Technicolor," Ph.D. Thesis, larXiv:0704.3999l [hep-ph]; 

K. Jan sen, J. Kuti and C. Liu , Phys. Lett. B 309, 119 

(1993) [arXiv:hep-lat /9305003| . 
[10] N. Nakanishi, Phys. Rev. D 3, 3235 (1971); idem, Phys. 

Rev. D 3, 811 (1971); idem, Phys. Rev. D 5, 1968 (1972); 

idem, Prog. Theor. Phys. Suppl. 51, 1 (1972). 
[11] S. Coleman, R. Jackiw, and H. D. Politzer, Phys. Rev. 

D 10 (1974) 2491, L. Abbott, J. Kang, and H. Schnitzer, 

Phys. Rev. D 13 (1976) 2212. 

[12] T. G. Rizzo, JHEP 0801 (2008) 042 |arXiv: 0712. 17911 

[hep-ph]], T. G . Rizzo, JHEP 0706 (2007) 070 

[arXiv:0704.3458l [hep-ph]] . 
[13] E. Alvarez, L. Da Rold, C. Schat and A. Szynkman, 

JHEP 0804 (2008) 026 |arXiv:0802. 10611 [hep-ph]]. 
[14] F. Wu and M. Z hong, Phys. Lett. B 659, 694 (2008) 

|arXiv:0705.3287l [hep-ph]] . 
[15] E. Tomboulis, Phys. Lett. B 70 (1977) 361, E. Tomboulis, 

Phys. Lett. B 97 (1980) 77. 



[16] K. Jansen, J. Kuti and C. Liu, Nucl. Phys. Proc. 
Suppl. 30, 681 (1993); idem, P hys. Lett. B 309, 
127 (1993) [arXiv:hep-lat/9305004| ; C. Liu, K. Jansen 
and J. Kuti, Nucl. Phy s. Proc. Suppl. 42, 630 (1995) 
[arXiv:hep-lat / 9412034] . 

[17] G. 't Hooft, Nucl. Phys. B 75 (1974) 461. 

[18] G. 't Hooft and M. J. G. Veltman, NATO Adv. Study 
Inst. Ser. B Phys. 4 (1974) 177. 

[19] Higher derivative gravity, a theory in which terms 
quadratic in the curvature tensor are added to the usual 
Einstein-Hilbert action, is known to be renormalizabc 
and asymptotically free. Partly for the reasons that mo- 
tivate this paper, it is not known whether this theory can 
be consistently formulated to all orders in perturbation 
theory. 

[20] Liu does check explicitly that the optical theorem holds 
for the scattering of Goldstone bosons in the sponta- 
neously broken version of his theory. 

[21] In fact we will need all inverse masses and momenta to 
be small compared to all the separations involved in the 
problem, including eventually the distance between col- 
lision and decay points, z' — zq. By restoring Planck's 
constant, |?7i| 3> h/m, we see that a simpler, equiva- 
lent assumption is that we analyze the kinematics in the 
semiclassical limit. We emphasize, however, that the ex- 
act scattering amplitude is used in this section, that is, 
no approximation for the scattering need be made in dis- 
cussing time dependence in a scattering experiment. 

[22] As explained above, this is justified by considering h — » 0. 

[23] Only the exponential has a fast oscillation with the vari- 
able q as h — > 0; the remaining factor does contain an 
exponential involving 1/h, namely exp(— im(j\ + T2)/h), 
however this is g-independent. 

[24] The flavor structure is suppressed here. 

[25] Except for one loop self-energy graphs, which give a triv- 
ial log divergence mass shift. 



